1/1 


AD- A 1 3 1  912  INTERACTION  OF  A  DISLOCATION  WITH  A  CRACK(U)  PRINCETON 
UN IV  Nij  DEPT  OF  CIVIL  ENGINEERING  A  C  ERINGEN  AUG  83 
RR-83-SM-9  N00014-76C-0240 


UNCLASSIFIED 


F/G  12/1 


NL 


OTIC  FILE  COPY 


Interaction  of  a  Dislocation  with  a  Crack 


A.  Cemal  Eringen 
PRINCETON  UNIVERSITY 


Technical  Report  No.  59 
Civil  Engng.  Res.  Rep.  No.  83-SM-9 


AUG  3  0  1983 


Research  Sponsored  by  the 
OFFICE  OF  NAVAL  RESEARCH 
under 

Contract  N00014-76-C-0240  Mod  4 
Task  No.  NR  064-410 


Aug.  1983 

Approved  for  public  release: 
distribution  unlimited 

Reproduction  in  whole  or  in  part  is  permitted 
for  any  purpose  of  the  United  States  Government 


08  26  ©64 


Interaction  of  a  Dislocation  with  a  Crack 

A.  Cemal  Eringen 
Princeton  University 
Princeton,  NJ  08544 

ABSTRACT 

A  solution  is  given  of  the  field  equations  of  nonlocal 
elasticity  for  a  line  crack  interacting  with  a  screw  dislocation  in 
an  elastic  plane  under  anti-plane  shear  loading.  Displacement  and 
stress  fields  are  determined  throughout  the  core  region  and  beyond. 

In  the  case  when  the  dislocation  is  absent,  the  circumferential  stress 
is  shown  to  vanish  at  the  crack  tip,  increasing  to  a  maximum  along  the 
crack  line  afterwards  decreasing  to  its  classical  value  at  large  dis¬ 
tances  from  the  crack  tip.  This  is  in  contradiction  with  the  classical 
elasticity  solutions  which  predicts  stress  singularity  at  the  crack 
tip  and  it  is  in  accordance  with  the  physical  condition  that  the  crack 
tip  surface  must  be  free  of  surface  tractions.  The  presence  of  the 
dislocation  alters  the  stress  distribution  considerably  when  it  is 
close  to  the  crack  tip.  The  stress  distributions,  in  the  core  region, 
are  displayed.  A  fracture  criterion  based  on  the  maximum  stress  is 
established  and  used  to  determine  the  theoretical  strengths  of  pure 
crystals  that  contain  a  line  crack.  Results  are  in  good  agreement  with 
those  based  on  the  atomic  theories  and  experiments. 
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ABSTRACT  ( Coni  mo r  mi  rarataa  lip  If  nacaaoary  and  Idonfly  A*  Bloch  nim>b»t' 

'A  solution  is  given  of  the  field  equations  of  nonlocal  elasti¬ 
city  for  a  line  crack  interacting  with  a  screw  dislocation  in 
an  elastic  plane  under  anti-plane  shear  loading.  Displacement 
and  stress  fields  are  determined  throughout  the  core  region  anc 
beyond.  In  the  case  when  the  dislocation  is  absent,  the  circun 
ferential  stress  is  shown  to  vanish  at  the  crack  tip,  increas- 
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(cont) 


ABSTRACT  (cont) 

to  its  classical  value  at  large  distances  from  the  crack 
tip.  This  is  in  contradiction  with  the  classical  elasticity 
solutions  which  predicts  stress  singularity  at  the  crack  tip 
and  it  is  in  accordance  with  the  physical  condition  that  the 
crack  tip  surface  must  be  free  of  surface  tractions.  The 
presence  of  the  dislocation  alters  the  stress  distribution 
considerably  when  it  is  close  to  the  crack  tip.  The  stress 
distributions,  in  the  core  region,  are  displayed.  A  frac¬ 
ture  criterion  based  on  the  maximum  stress  is  established 
and  used  to  determine  the  theoretical  strengths  of  pure 
crystals  that  contain  a  line  crack.  Results  are  in  good 
agreement  with  those  based  on  the  atomic  theories  and  ex¬ 
periments. 
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I.  INTRODUCTION 

It  is  well  known  that  the  classical  elasticity  solution  of  crack 
problems  fail  in  a  core  region  around  a  sharp  crack  tip,  since  they  pre¬ 
dict  stress  singularity  at  the  tip.  The  assessment  of  the  core  radius 
and  the  stress  field  within  the  core  is  a  problem  usually  discussed  within 
the  context  of  atomic  theories  of  lattices  (cf.  [1]),  even  at  that  its 
treatment  contains  various  assumptions  regarding  the  interatomic  arrange¬ 
ments  and  force  fields. 

Engineering  fracture  mechanics,  on  the  other  hand,  is  based  on  the 
2 

Griffith's  ideas  which  resort  to  other  concepts  (e.g.  energy,  J-integral, 
fracture  toughness).  To  be  sure,  there  exist  certain  erzatz  to  account 
for  the  effect  of  the  core  region  on  fracture  process  in  phenomenological 
ways.  These  are  useful  for  engineering  purposes,  however,  they  are  not 
based  on  a  fundamental  theory  nor  are  they  capable  predicting  the  state 
of  stress  in  the  core  region  whcih  is  fundamental  to  the  initiation  of 
fracture. 

In  several  previous  papers,  we  have  shown  that  nonlocal  elasticity 

solutions  of  Griffith  crack  problems  lead  to  finite  stress  at  the  crack 
3-5  6 

tip  .  In  fact,  an  exact  solution  obtained  for  the  screw  dislocation, 
indicates  that  the  stress  vanishes  at  the  tip  of  the  crack,  growing  to  a 
maximum  in  the  vicinity  of  the  crack  tip.  The  important  implications  of 
this  result  in  connection  with  the  initiation  of  fracture  is  the  motivation 
for  the  present  work. 

The  solution  obtained  here  for  the  Mode  III  (anti-plane  shear) 
problem  for  a  crack  interacting  with  a  screw  dislocation  indicates  that  the 
circumferential  stress  field  is  vanishingly  small  (zero  when  the  dislocation 


3 


is  absent)  at  the  crack  tip,  when  the  screw  is  located  far  away  from  the 
crack  tip.  When  the  dislcoation  is  near  the  crack  tip,  the  stress  field 
is  affected  apprciably  displaying  several  maxima  near  the  crack  tip.  By 
equating  the  maximum  stress  to  the  cohesive  yield  stress,  we  can  determine 
the  stress  intensity  factor  Kg  or  the  theoretical  yield  stress,  given 
Kg.  Calculated  ^Kg  values,  on  the  basis  of  the  present  theory,  are  in 
fair  agreement  with  those  determined  experimentally.  Theoretical  strengths 
are  also  estimated  by  means  of  the  dislocation  model.  Results  agree  with 
those  predicted  by  atomic  models. 

The  mathematical  model  of  approach  to  the  solution  of  this  problem 
is  new  and  possesses  potential  applications  in  other  areas. 


2.  BASIC  EQUATIONS 

In  several  previous  papers,  we  developed  a  theory  of  nonlocal 
elasticity,  cf.  [7,  8,  9],  For  homogeneous  and  isotropic  elastic  solids, 
linear  theory  is  expressed  by  the  set  of  equations 


(2.1) 


o  . 


(2.2) 

(2.3) 


(2.<) 


tkl(i.t) 

°kt(;’>t) 

• u <*'•*> 


•  J  o(|x'-x|,  T )  ok((x',t)  dv(x')  , 
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*  X  e 


rr(;'.t)  a 

3uk(x',t) 


k£ 
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+  2  «  •«<*’•»> 

3uf(x',t) 
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where  tki,  p  ,  f^  and  u£  are  respectively,  the  stress  tensor,  mass  density, 
body  force  density  and  the  displacement  vector,  A  and  y  are  the  Lam#  elastic 
constants  and  a  Is  the  "attenuation  function"  which  depends  on  the  distance 
1 x'-x |  and  a  parameter  t  which  denotes  the  ratio  of  the  internal  charac¬ 
teristic  length  a  to  the  external  characteristic  length  l  ,  i.e. 

(2.5)  t  *  eg  a/£ 

where  eg  is  a  constant  appropriate  to  each  material.  Characteristic 
lengths  may  be  selected  according  to  the  range  and  sensitivity  of  the 
physical  phenomena  to  be  investigated.  For  instance,  for  perfect  crystals, 
a  may  be  taken  as  the  lattice  parameter  and  £  as  the  half  crack  length. 

For  granular  materials,  a  may  be  considered  to  be  the  average  granular 
distance  and  for  fiber  composites,  the  fiber  distance  etc.  The  material 
constant,  eQ  may  be  determined  by  one  experiment. 

Equations  (2.1),  (2.3)  and  (2.4)  are  those  known  from  the  theory  of 
classical  elasticity,  but  Eq.  (2.2)  is  new,  replacing  Hooke's  law.  Accord¬ 
ing  to  Eq.  (2.2),  the  stress  at  a  point  x  depends  on  strains  at  all  points 
x'  of  the  body.  The  attenuation  function  a  determines  the  degree  of 
influence  with  the  distance.  From  the  physical  nature  of  solids,  it  is 
clear  that  the  influence  of  strains  at  x*  ,  on  the  stress  at  x  ,  decreases 
with  the  distance  Ix'-xl  .  Thus,  a(|x'-x!)  must  acquire  its  maximum  at 
x'«x.  Moreover,  when  a -*■(),  a  must  become  a  Dirac  delta  measure  so  that 
nonlocal  theory  shall  revert  to  classical  elasticity  theory.  By  matching  the 
phonon  dispersion  curves  with  those  resulting  from  nonlocal  theory,  we  have 
determined  a  for  various  cases  (cf.  [5],  [83*  P03)* 
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i 

By  discretizing  Eq.  (2.2),  it  can  be  shown  that  equations  of  nonlocal 

1  1  ^ 
elasticity  revert  to  those  of  atomic  lattice  dynamics.  Thus,  it  is  clear  j 

that  nonlocal  theory  is  a  suitable  model  for  the  treatment  of  physical  pheno¬ 
mena  with  characteristic  lengths  in  the  range  from  the  molecular  or  atomic 
dimensions  to  macroscopic  sizes. 

For  a  two-dimensional  perfect  lattice,  the  dispersion  curves  are 

matched  in  the  entire  Brillouin  zone  to  within  an  error  less  than  1%« 

10 

with  the  attenuation  function 

(2.6)  a(  |x| ,  t  )  =  (2ti£2t2)”1  K^^yx/Z-r) 

where  KQ  is  the  modified  Bessel's  function.  We  note  that  Eq.  (2.6)  is 

2  2  2 

Green's  function  for  the  operator  L  =  (1  -  £  t  V‘  ),  i.e. 


(2.7)  (1  -  £2t2V2)o  =  6( I x ' -x | ) 

In  fact,  it  is  possible  to  employ  other  linear  operators  to  characterize  the 
nature  of  nonlocal  attractions  of  material  points  in  solids.  This  apparent 
non-uniqueness  of  a  may  be  considered  to  be  a  defect  of  the  theory.  On 
the  contrary,  for  imperfect  and  amorphous  solids,  this  may  provide  a  de¬ 
sirable  flexibility.  Ultimately,  however,  a  should  be  determined  from 
experimental  and/or  statistical  mechanical  considerations.  For  perfect 
crystals,  Eq.  (2.6)  leads  to  excellent  agreements  with  the  dispersion  curves 
based  on  atomic  lattice  theory. 


■rap 
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(2.8)  (1  -  (W)  tu  ■  0|(i 

Divergence  of  Eq.  (2.8),  upon  using  (2.1)  and  (2.3),  leads  to 


(2.9)  (X+u)  uk^k£  +  yu1)k|(  +  (1  -  £2i2V2)(pf£  -  pu£)  *  0 

For  the  static  case  and  vanishing  body  forces,  Eq.  (2.9)  is  non  other  than 
Navier's  equation  of  classical  elasticity.  Note,  however,  that  the  stress 
tensor  is  not  okJ.  but  tk£  and  it  requires  that  we  solve  Eq.  (2.8)  to 
determine  tk£  . 

For  plane,  harmonic,  SH-waves,  Eq.  (2.9)  gives  the  frequency 

(2.10)  u  *  (y/p)Ji  k[l  +  e2k2a2r* 

where  k  is  the  amplitude  of  the  wave  vector.  By  equating  w  given  by 
Eq.  (2.10)  to  that  predicted  by  the  Bom-K3rmSn  model  of  lattice  dynamics, 
at  the  end  of  the  Brillouin  zone  (ka*Tt),  we  find  that 

(2.11)  eQ  -  (tt2  -  4)*/2ir  «  0.39  . 

The  dispersion  curve  based  on  Eq.  (2.10)  and  that  of  the  Bom-KSrmSn  model 
are  compared  in  Fig.  1.  We  see  that  the  matching  is  very  good.  The  maximum 
error  is  less  than  6*. 


0 
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The  dispersion  curves  of  the  Bom-KirmSn  model  is  a  good  approximation 
for  some  fee  and  bee  metals  (e.g.  A1  and  cu).  While  the  Brillouin 
zone  may  vary  In  different  directions  of  slips  in  crystals,  we  believe  that 
Eq.  (2.11)  is  a  reasonable  value  for  eQ  when  the  material  Is  considered 

to  be  isotropic. 


3.  CLASSICAL  STRESS  FIELDS 

A  homogeneous,  isotropic  elastic  solid  of  infinite  extent  contains 
a  crack  located  at  -  c  £  x,  <_  c  ,  x3  =  0  ,  -®  <x3 < ®  where  xk  are  the 
rectangular  coordinates.  Fig.  2.  We  suppose  that  there  exists  a  disloca¬ 
tion  which  lies  parallel  to  the  x3-axis  and  which  intersects  the  plane  x3=0 
at  the  point  S{x-|  =  £  ,  x2  =  n)  .  The  solid  is  subject  to  a  constant  anti¬ 
plane  shear  at  x^ =  +  ®  .  The  classical  elasticity  solution  of  this  problem 
12 

was  given  by  Louat  .  However,  here  we  derive  the  solution  of  this  problem 
in  the  form  better  suited  for  our  purpose,  eliminating  possible  misprints, 
difficulties  in  notations  and  in  taking  various  limits. 

Since  the  state  of  the  body  is  the  same  at  all  planes,  x3= const., 
the  problem  is  two-dimensional  and  we  need  to  treat  the  plane  problem  in  the 
plane  x3=0  with  a  1  ine  crack  located  at  X£  s  0  ,  |x^|£C. 

The  classical  stress  field  at  any  point  Pfx^.x^)  may  be  expressed 


conveniently  in  the  form 
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(3.1) 


where  z  =  x1  -  i  and  f(t)  is  the  distribution  function  which  is  the 
solution  of  the  equation  of  equilibrium  of  the  forces  acting  on  the  crack 
surface: 


*c 

(3.2)  A  f  dt  =  o^x)  +  Oq  ,  A  =  uXq/2tt 

-c 

Here,  the  integral  denotes  a  Cauchy  principal  value,  y  is  the  shear  modul 
Aq  is  the  displacement  vector  of  a  unit  positive  dislocation  and  oQ  and  oc 
are  the  stress  fields  at  the  crack  surface  due  to  the  applied  load  and  the 
dislocation,  respecti vely. 

The  solution  of  the  integral  equation  (3.2)  is  well-known,  Tricomi13 


(3.3)  f(x)  =  - 


rV 


:2-x2 


/ 


-c 


[Of 


od(t)] 


dt 

t-x 


Here,  Q  is  a  constant  to  be  determined  from  the  condition  that 


A 


(3.4) 


n 


where  nXQ  is  the  total  dislocation  content  of  the  distribution  f(x)  . 
The  stress  o^t)  is  given  by 
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(3.5' 


°d(t) 


—  Vb(ty,) 
2tt[U-02  +  n3 


Substituting  this  into  Eq.  (3.3),  we  can  carry  out  integrations  to 

obtain 


f(x)  = 


(3.6) 


*?; 


I°0X 


a  r' 

TtX,  L 


~2  T 


2(;-x) 


2(;-x) 


1]  +  q} 


where  ;  =  £  +  in  ,  ;  =  £  -  in  .  Using  this  in  Eq.  (3.4),  we  will  have 

Q  =  n/~  .  Carrying  f(x)  into  Eq.  (3.1)  after  some  tedious  integrations, 
we  obtain 


(3.7) 


c23  "  1  °13  =  °0(  / 


J=-  i)  -f.tr1-  <> 

=2  _2  ZX0  Z-C 


z  -c 


Z-C 


1  b  A 

♦  ~r  0  -  4-=-)]  +  (i+  n)  ~— 


z  -c 


z  -c 


The  forces  acting  on  the  dislocation  at  (£,ri)»  due  to  the  crack, 
are  given  by 

(3.8)  F-|  -  b  022*  ^  2  ~  b  ^  ( x^  -  t,  *  X£  ~  n ) 


For  our  own  purpose  later,  we  need  the  stress  field  when  the  dislocation  is 
located  along  the  x-j-axis  and  the  surface  of  the  crack  is  free  of  tractions. 
To  this  end,  we  set  n = 0  and  add 
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(3.9) 


°o  +  0d(xi)  “  0o+>,rrr^T 


to  the  right-hand  side  of  Eq.  (3.7).  Hence, 


(3.10) 


°23  •  '  °13  +  °0  +  °d(xl) 


— - —  [oD  z  +  A(^  +  n) 

^ r; j  0  xo 


+  bA  52] 

A0  z-c 


gives  the  classical  stress  field  at  any  point  outside  of  the  crack  when 
the  body  is  loaded  at  X£  =  +  *  with  a  constant  shear  0^3  =  —  °D  '  When 
the  crack  contains  no  dislocations,  then  we  have  n  =  0. 

Two  special  cases  are  important: 

(i)  Cr2?<  srJ  0^  =  0.  In  this  case,  the  classical  stress  field 

is  given  by 

(3.11)  o  «  -*&- 

2tt  z 


where  we  also  set  £  =  0  placing  the  dislocation  to  the  origin  of  coordinates. 


(ii)  Kc  Dislocations.  In  this  case,  A=0  and  we  have 


(3.12) 


v 


Both  of  these  results  are  well-known  in  the  literature. 
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4.  NONLOCAL  STRESS  FIELDS 

To  determine  the  nonlocal  stress  fields,  we  must  obtain  the  solution 
of 

(4.1)  (l-t2£2V2)t  =  o 
subject  to  some  boundary  conditions.  Here 

(4.2)  t  =  t23  -  i  t^  ,  o  =  o23  -  i  o13  +  oQ  +  od(x-|) 

2 

Since  V  o  =  0,  t  =  o  is  a  particular  solution  of  Eq.  (4.1).  The  comple¬ 
mentary  solution  of  (4.1),  vanishing  at  infinity  and  having  proper  symmetry 
regulations  with  respect  to  (xj,  +  X2)  »  1S  of  form 

(4.3)  tc  =  Kv(r/iT)(Av  e1v6  +  Bv  e'iv9) 

where  A  ,  B  and  v  are  constants,  K  (p)  is  the  modified  Bessel's 

v  v  v 

function  and  (r,6)  are  the  plane  polar  coordinates. 

The  boundary  condition  on  the  crack  surface  requires  that  t23=0. 
Taking  the  origin  ra0  of  the  coordinates  at  the  right-hand  crack  tip 
and  writing  r =  r^  ,  6 = 0^  ,  in  (4.3)  we  see  that  to  fulfill  this  condi¬ 

tion,  we  must  have  v=l/2  ,  since  all  other  solutions  lead  to  displacement 
singularities  at  r-j  =  0  . 

Classical  stress  field  o  possesses  singularity  at  the  screw 
dislocation  =  £  ,  x2  *  0  .  The  surface  traction,  trz  on  the  edge 
surface  of  the  dislocation  is  required  to  vanish,  according  to  the  boundary 
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condition.  To  fulfill  this  condition  we  take  v  =  l  and  move  the  origin 
of  coordinates  to  x^  *£,  x2 = 0  .  This  may  be  expressed  by  writing 
r  =  rd  ,  6  =  63  and 

ie,  ie, 

(4.4)  rd  e  =  r]  e  '  xo  • 

Hence,  the  general  solution  of  (4.1)  appropriate  to  our  problem  is  of  the 
form 

,  -r,/t£  ie,/2  -ie-,/2 

(4.5)  t  =  (mZ/Zr^P  e  (C1  e  '  +  C2  e  ) 

ie,  -ie, 

+  K1(rd/t£)(C3  e  J  +  C4  e  J)  +  o  . 


To  determine  C&  ,  we  calculate  stress  components  in  polar  coordinates 
(r,,e,)  : 

(4'6)  hz  '  1  Vz  *  (t23  ‘  1  »13>  e  1  • 


We  imagine  the  crack  tip  as  a  limit  of  a  small  circular  arc  with  radius 
r1  = e  approaching  zero.  For  small  c  ,  we  have  approximately 


(4.7) 


z 

Z1 


c  +  z. 


i©-. 


e  e 


-  c  +  z2  =  C  +  23  =  c  +  xQ  +  z3  , 


Using  these  in  Eq.  (3.10),  we  will  have 


13 


Next,  we  calculate  the  stress  field  at  the  location  rrf  =  0  of  the  screw 
dislocation.  As  rd+0  we  have  r^  =  Xq  and 

*In  fact,  this  condition  is  satisfied  exactly  along  the  crack  line 
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v  -xJxsl  ie,  -ie, 

(4.12)  tze  -  i  t2r  -  (ttt£/2x0P  e  u  C1  +  K1(rd/T£)(C3  e  J  +  C4e  J 

+  [x0(x0+2c)]'V0(c+x0)  +  ^-o+-r) 

+  2iHX0(x0  +  2c^  e  ^  * 

d 

Again  tzr  must  vanish  as  rd -*•().  This  implies  that  C4  s  0  and 

(4.13)  C3  =  -yb/27t  t£  . 

The  general  solution  is  now  complete. 

,  -r./ii,  ie,/2  ie, 

(4.14)  t  =  (7TT£/2r1)^  e  e  +  K-j(rd/x£)  C3  e  J  +  o 


where  C-|  and  C3  are  given  by  (4.11)  and  (4.13).  In  polar  coordinates, 
we  have 


(4.15) 


,  ” r.  /t £  ~ i  6 1  /2  i  C  0 "a- 8 1 ) 

t02-i  trz  =  (irT£/2r1)«  e  e  '  +  K^/xt)  C3  e  J  1 

.  /  1  (e2"el ^/2/r  -ie  ,ub  M  .  nX0> 

+  ( rl  r2 )  5  e  j[oQ  re  0  +  — ) 

+  $  (ri  e"61-  V1  [«o(«o  * 


Special  cases  mentioned  in  Section  3  can  be  obtained  in  a  similar  fashion. 
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(i )  No  Crack 


(4.16) 


where  we  have  taken  the  origin  of  the  polar  coordinates  at  the  dislocation, 
so  that  p  =  r^/t i  . 

( i i )  No  Dislocation 


(4.17) 


tez " 1  Vz 


o0(c/2ri)*  [(2r2/cr2)^  e 


H-e  +  y-) 


-ie1/2 

e  ]  e 


5.  FRACTURE 

Here,  we  discuss  the  onset  of  fracture  and  determine  the 
theoretical  stresses  for  the  two  special  cases. 

(i )  No  Crack 

According  to  Eq.  (4.16),  we  have  trz *  0  and 
(5.1)  Te(p)  •  t,,  ■  p-'[1-pK,(p)] 


The  maximum  of  T0  occurs  at  p  « 1 . 1  and  is  given  by 


T 


©max 


1.1 


(5.2) 


0.3993  ; 
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It  is  natural  to  assune  that  when  ta_mmw  becomes  equal  to  the  theoretical 

ez  max 

stress  ty  ,  the  crystal  will  rupture.  Thus, 

(5.3)  V„  =  0-3993  ^ 

If  we  write  h * e0a/0.3993,  this  agrees  with  the  estimate  of  Frenkel  based 
on  an  atomic  model  (cf.  Kelly^,  p.  12).  For  aluminum  (fee),  b  =  a/>/?  and 
for  iron  (bcc),  b=  /3/2,  so  that  Eq.  (5.3)  gives 


0.12 

{AT :  [111]  <1T0  >} 

0.14 

(Fe:  [lTo] < 111>) 

These  are  close  to  the  theoretical  results  ty/p  =  0.11  based  on  atomic 
models. 

It  is  interesting  to  note  that  t02 *  0  at  the  center  of  dislocation 
and  it  rises  to  a  maximum  at  p*  1.1  ,  thereafter  decreasing  to  zero  with 
p  .  Significant  consequences  of  the  present  predictions  as  contrasted  to 
the  classical  results  are: 

(a)  The  stress  at  the  center  of  the  core  is  not  infinite,  but  zero. 

(b)  Fracture  begins  at  p  *  p£  not  at  the  center  of  the  core. 

(c)  There  is  a  low  stress  region,  0<p<pc  within  the  core. 

(11)  No  Dislocations 

From  (4.17),  it  is  clear  that  t02  acquires  its  maximum  along  the  crack 
line  e«  e1  «  e  2  *0,  near  the  crack  tip.  The  circumferential  stress  along 
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the  crack  line  >  0  is  expressed  by 

(5.4)  tz6^°0  “  (2yp)”^  [(1  +  yp) ( 1  +  -y)  ^  -  e  p] 
where 

(5.5)  p  =  r^ /eQa  ,  y  c  eQa/c 

t  q  vanishes  at  the  crack  tip  p*0  and  has  a  maximum  at  pspc  which  is 
the  root  of 

(5.6)  e"p  (1  +  2p)  -  (1  +  ^)“3/2 

Since  Y«1  (>  <  10"6),  we  see  that  the  root  of  (5.6)  is  independent  of 

c  and  is  given  by 

(5.7)  Pc  =  1.2565 
and  the  maximum  stress  is  given  by 


(5'8>  ‘ze  max  *  °0(e0a/c,'!S  ^  • 

We  also  observe  that  as  p -*•«>,  (5.4)  gives  t  g  *  °q  ’  as  **  should. 

In  the  classical  tradition,  if  we  write  Kjjj  *  /ire  Og  ,  then  (5.4) 


may  be  expressed  as 


k 
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(5.9)  T0(p)  =  A  t20/oo  =  (it  eQa)^  tze/Km 

=  (2p)‘*  [1+yp)(1  +  ^)‘%  -  e’p] 


This  is  plotted  against  p  ,  in  Fig.  3  in  the  vicinity  of  the  crack  tip. 
The  classical  (local)  stress  is  also  indicated  on  this  figure  by  dashed 
lines.  From  this  figure,  it  is  clear  that  the  classical  stress  field 
deviates  considerably  from  the  nonlocal  stress  field  in  the  region 
0  £  p<  5.  In  fact,  it  diverges  at  the  crack  tip. 

A  perfect  crystal  which  contains  a  crack,  but  no  dislocation,  will 
not  rupture  before  the  maximum  stress  reaches  the  value  of  the  cohesive 
stress  (theoretical  stress)  that  holds  the  atomic  bonds  of  the  lattice. 
Thus,  the  entire  crystal  is  in  the  elastic  state  of  equilibrium  when 


(5.10) 


tz6max  <  *y 


The  failure  begins  when  tz0max  =  t^  ,  i.e. ,  when 


(5.11) 


Kc/ty  =  (*  e0a)*  (</2p c  *  — ~ ) 


ftp ~ 


3.9278 


where  Kc  =  /ire  Cqc  is  the  critical  fracture  toughness. 

Using  Eq.  (5.11)  and  eQ*0.39,  we  calculate  a  few  Kg-values  which 
are  listed  in  Table  1  (last  column)  along  with  classical  Kc's  based  on 
Kc  *  4pys)^,  where  is  the  surface  energy.  Experimental  observations 
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15,16 

of  Ohr  and  Chang  are  also  listed  in  this  Table.  Classical  estimates 

are  expected  to  be  inaccurate  considering  the  fact  that  even  with  the  best 
present-day  techniques  available,  they  could  not  be  measured  to  an  accuracy 
better  than  a  factor  of  two.  Moreover,  classical  formula  assumes  no  de¬ 
fects  (i.e. ,  no  crack  and  dislocation),  therefore  it  is  expected  to  give 
higher  Kc*values.  On  the  other  hand,  experimental  measurements  of  Ohr 
and  Chang  required  the  measurement  of  the  length  of  the  plastic  zone  among 
other  constants.  This  already  implies  the  existence  of  dislocations  so 
that  we  expect  some  deviation  from  the  perfect  crystal  containing  no  dis¬ 
location  but  a  single  crack.  In  Section  6,  we  examine  the  general  case 
when  the  solid  contains  a  crack  and  a  dislocation. 


6.  DISLOCATION  AND  CRACK 

Along  the  crack  line  x^  >  c,  x^  =  0  , 

Eq.  (4.15)  gives  t  =  0  and  with  n  =  0 

(fi.l)  tez  -  tc  ♦  tdc 

where 

(6.2)  tc  ST  /  °o  =  T]  =  (2o)_?J  [1  +  yp) (1  +  -  e"c], 

tdc  H  /o0B  =  T2  =  (2 c)'h  {[O^r*4  -  [1  -(1+  4-  )S': 

-  sgn(o-x  )(2p/y  )  *  K^jo-i  |)  +  (1+  f-  )'H  (n4  )\f-  -l)'1 

°  '  0  0 


(6.3) 
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in  which 


(6.4) 


xo  = 


w 


when  the  dislocation  is  absent  we  have  tdc  =  0  so  that  tdc  is  the 
shear  stress  arising  from  the  interaction  of  the  dislocation  with  the 
crack.  At  the  crack  tip  p  =  0  and  we  have 


(6.5) 


Hz' 


V 


=  -  K. 

y  l 


(*n> 


This  shows  that  for  positive  dislocation  (  b>0  )  the  shear  stress  is 
positive  and  therefore  the  crack  tip  will  tend  to  close  up  for  b>0. 

For  p  <  0  the  opposite  will  occur.  However,  the  stress  given  by  (6.5) 
due  to  the  dislocation  is  very  small  for  large  Xq  and  it  becomes  large 
when  the  dislocation  is  very  close  to  the  crack  tip.  The  stress  field 
T-j  ( .  )  =  T^_ ( ; )  so  that  Fig.  3  represents  T^(p)  .  Fig.  4  displays 
graphs  of  T^;)  for  several  values  of  xQ  =  3.1,  4  1,  6.1  and  10.1, 

O 

keeping  •  =  10  fixed.  These  graphs  show  that  T^  possesses  a  minimum 
and  two  maxima.  The  crack  tip  is  not  stress-free.  The  maximum 
of  T2  occur  close  to  the  dislocation,  For  example,  in  the  case 
x  =  3.1,  the  maximum  is  at  xc  =  4.4  and  in  the  case  xQ  =  10.1  it  is 
at  x  =  1 1 .2. 

To  obtain  an  idea  on  the  combined  effect  I  have  selected 
-4 

c  =  10  and  plotted  the  ratio  of  the  combined  stress  to  00(tez/'V 

xo  .  From  these  graphs  it  is  clear  that  when 


in  Fig.  5  for  various 
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the  dislocation  is  located  close  to  the  maximum  of  the  combined 

effect  is  large.  For  example,  while  T,  =  0.451  for  the  case  of 

lmax 

crack  alone,  for  the  combined  effect  we  have  T  =  0.63  so  that 

max 

the  ratio  of  the  two  K  - values  is  given  by 


(6.6) 


Vot  '  Kgi 


0.71 


This  implies  that  a  dislocation  located  at  a  distance  approximately  one 
or  two  lattice  parameters  away  from  the  crack  tip  reduces  the  fracture 
toughness  by  about  30°i.  Hence  the  theoretical  values  of  listed 

in  Table  1  will  be  reduced  about  30«  bringing  the  numbers  on  the  last 
column  closer  to  those  listed  in  the  adjacent  column  marked  experiments. 
These  results  however  must  still  be  considered  only  as  indications  for 
the  trend.  A  more  realistic  physical  picture  requires  the  presence  of 
large  numbers  of  dislocations  distributed  over  a  few  microns  or  so 
distance,  away  from  the  crack  tip.  Consequently,  to  obtain  a  close 
approximation  to  experimental  observations  of  Chang  and  Ohr  we  need  to 
consider  a  distribution  of  dislocation  in  a  region  near  the  crack  tip. 
Such  a  consideration  will  require  a  separate  study  of  dislocation  pile  up 
which  is  left  to  a  future  study. 
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and  Structures 
OUSDR4E,  The  Pentagon 
Soon  3D1089 

Washington,  D.C.  20301 

Chief,  Airframe  and  Equipment  Branch 
FS-120 

Office  of  Flight  Standards 
Federal  Aviation  Agency 
Washington,  D.C.  20553 

National  Academy  of  Sciences 
National  Research  Council 
Ship  Hull  Research  Committee 
2101  Constitution  Avenue 
Washington,  D.C.  20418 
Attn:  Mr.  A.  R.  Lytle 

National  Science  Foundation 
Engineering  Mechanics  Section 
Division  of  Engineering 
Washington,  D.C.  20550 

Picatinny  Arsenal 

Plastics  Technical  Evaluation  Center 
Attn:  Technical  Information  Section 
Dover,  New  Jersey  07801 

Maritime  Administration 
Office  of  Maritime  Technology 
14th  and  Constitution  Avenue,  NW. 
Washington,  D.C.  20230 

PART  2  -  Contractors  and  Other  Technical 
Collaborators 

Universities 

Dr.  J.  Tinsley  Oden 
University  of  Texas  at  Austin 
345  Engineering  Science  Building 
Austin,  Texas  78712 

Professor  Julius  Miklovitz 
California  Institute  of  Technology 
Division  of  Engineering 
and  Applied  Sciences 
Pasadena,  California  91109 


Universities  (Con't) 

Dr.  Harold  Liebovitz,  Dean 
School  of  Engineering  and 
Applied  Science 
George  Washington  University 
Washington,  D.C.  20052 

Professor  Eli  Sternberg 
California  Institute  of  Technology 
Division  of  Engineering  and 
Applied  Sciences 
Pasadena,  California  91109 

Professor  Paul  M.  Naghdi 
University  of  California 
Department  of  Mechanical  Engineering 
Berkeley,  California  94720 

Professor  A.  J.  Durelli 
Oakland  University 
School  of  Engineering 
Rochester,  Missouri  48063 

Professor  F.  L.  DiMaggio 
Columbia  University 
Department  of  Civil  Engineering 
New  York,  New  York  10027 

Professor  Norman  Jones 

The  University  of  Liverpool 

Department  of  Mechanical  Engineering 

P.  0.  Box  147 

Brown low  Hill 

Liverpool  L69  3BX 

England 

Professor  E.  J.  Skudrzyk 
Pennsylvania  State  University 
Applied  Research  Laboratory 
Department  of  Physics 
State  College,  Pennsylvania  16801 

Professor  J.  Klosner 
Polytechnic  Institute  of  New  York 
Department  of  Mechanical  and 
Aerospace  Engineering 
333  Jay  Street 
Brooklyn,  New  York  11201 

Professor  R.  A.  Schapery 
Texas  ASM  University 
Department  of  Civil  Engineering 
College  Stetion,  Texes  77843 
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Universities  (Con't.) 

Professor  Walter  D.  Pilkey 
University  of  Virginia 
Research  Laboratories  for  the 
Engineering  Sciences  and 
Applied  Sciences 
Charlottesville,  Virginia  22901 

Professor  K.  D.  Willaert 
Clarkson  College  of  Technology 
Department  of  Mechanical  Engineering 
Potsdam,  Mew  York  13676 

Or.  Walter  E.  Haisler 
Texas  ASM  University 
Aerospace  Engineering  Department 
College  Station,  Texas  77843 

Dr.  Bussein  A.  Kamel 
University  of  Arizona 
Department  of  Aerospace  and 
Mechanical  Engineering 
Tucson,  Arizona  85721 

Dr.  S.  J.  Fenves 
Camegie-Mellon  University 
Department  of  Civil  Engineering 
Schenley  Park 

Pittsburgh,  Pennsylvania  15213 

Dr.  Ronald  L.  Huston 
Department  of  Engineering  Analysis 
University  of  Cincinnati 
Cincinnati,  Ohio  45221 

Professor  G.  C.  M.  Sih 
Lehigh  University 
Institute  of  Fracture  and 
Solid  Mechanics 
Bethlehem,  Pennsylvania  18015 

Professor  Albert  S.  Kobayashi 
University  of  Washington 
Department  of  Mechanical  Engineering 
Seattle,  Waahington  981Q5 

Professor  Daniel  Frederick 
Virginia  Polytechnic  Institute  and 
State  University 

Department  of  Engineering  Mechanics 
Blacksburg,  Virginia  24061 


Universities  (Con't) 

Professor  A.  C.  Eringen 
Princeton  University 
Department  of  Aerospace  and 
Mechanical  Sciences 
Princeton,  Mew  Jersey  08540 

Professor  E.  H.  Lee 

Stanford  University 

Division  of  Engineering  Mechanics 

Stanford,  California  94305 

Professor  Albert  I.  King 
Wayne  State  University 
Biomechanics  Research  Center 
Detroit,  Michigan  48202 

Dr.  V.  R.  Hodgson 
Wayne  State  University 
School  of  Medicine 
Detroit,  Michigan  48202 

Dean  B.  A.  Boley 
Northwestern  University 
Department  of  Civil  Engineering 
Evanston,  Illinois  60201 

Professor  P.  G.  Hodge,  Jr. 
University  of  Minnesota 
Department  of  Aerospace  Engineering 
and  Mechanics 

Minneapolis,  Minnesota  55455 

Dr.  D.  C.  Drucker 
University  of  Illinois 
Dean  of  Engineering 
Urbane,  Illinois  61801 

Professor  M.  M.  Newmark 
University  of  Illinois 
Department  of  Civil  Engineering 
Urbane,  Illinois  61803 

Professor  E.  Reissner 
University  of  California,  San  Diego 
Department  of  Applied  Mechanics 
La  Jolla,  California  92037 

Professor  William  A.  Mash 
University  of  Massachusetts 
Department  of  Mechanics  and 
Aerospace  Engineering 
Amherst,  Massachusetts  01002 
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Universit ies  (Con't) 

Professor  G.  Herrmann 
Stanford  University 
Department  of  Applied  Mechanics 
Stanford,  California  94305 

Professor  J.  D.  Achenbach 
Northwest  University 
Department  of  Civil  Engineering 
Evanston,  Illinois  60201 

Professor  S.  B.  Dong 
University  of  California 
Department  of  Mechanics 
Los  Angeles,  California  90024 

Professor  Burt  Paul 
University  of  Pennsylvania 
Towne  School  of  Civil  and 
Mechanical  Engineering 
Philadelphia,  Pennsylvania  19104 

Professor  H.  W.  Liu 
Syracuse  University 
Department  of  Chemical  Engineering 
and  Metallurgy 
Syracuse,  New  York  13210 

Professor  S.  Bodner 
Technion  RAD  Foundation 
Haifa,  Israel 

Professor  Werner  Goldsmith 
University  of  California 
Department  of  Mechanical  Engineering 
Berkeley,  California  94720 

Professor  R.  S.  Rivlin 
Lehigh  University 
Center  for  the  Applies  ion 
of  Mathemat ics 

Bethlehem,  Pennsylvania  18015 

Professor  F.  A.  Cozzarelli 
State  University  of  New  York  at 
Buffalo 

Division  of  Interdisciplinary  Studies 
Karr  Parker  Engineering  Building 
Chemistry  Road 
Buffalo,  New  York  14214 


Universit ies  (Con't) 

Professor  Joseph  L.  Rose 
Drexel  University 

Department  of  Mechanical  Engineering 
and  Mechanics 

Philadelphia,  Pennsylvania  19104 

Professor  B.  K.  Donaldson 
University  of  Maryland 
Aerospace  Engineering  Department 
College  Park,  Maryland  20742 

Professor  Joseph  A.  Clark 
Catholic  University  of  America 
Department  of  Mechanical  Engineering 
Washington,  D.C.  20064 

Dr.  Samuel  B.  Batdorf 
University  of  California 
School  of  Engineering 
and  Applied  Science 
Los  Angeles,  California  90024 

Professor  Isaac  Fried 
Boston  Universitv 
Department  of  Mathematics 
Boston,  Massachusetts  02215 

Professor  E.  Krempl 
Rensselaer  Polytechnic  Institute 
Division  of  Engineering 
Engineering  Mechanics 
Troy,  New  York  12181 

Dr.  Jack  R.  Vinson 
University  of  Delaware 
Department  of  Mechanical  and  Aerospace 
Engineering  and  the  Center  for 
Composite  Materials 
Newark,  Delaware  19711 

Dr.  J.  Duffy 
Brown  University 
Division  of  Engineering 
Providence,  Rhode  Island  02912 

Dr.  J.  L.  Swedlow 
Camegie-Mellon  University 
Department  of  Mechanical  Engineering 
Pittsburgh,  Pennsylvania  15213 
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Universities  (Con't)  Pniversit ies  (Con't) 


Dr.  V.  K.  Versdan 

Ohio  State  University  Research  Foundation 
Department  of  Engineering  Mechanics 
Columbus,  Ohio  43210 

Dr.  Z.  Hash in 
University  of  Pennsylvania 
Department  of  Metallurgy  and 
Materials  Science 
College  of  Engineering  and 
Applied  Science 

Philadelphia,  Pennsylvania  19104 

Dr.  Jackson  C.  S.  Yang 
University  of  Maryland 
Department  of  Mechanical  Engineering 
College  Park,  Maryland  20742 

Professor  T.  Y.  Chang 
University  of  Akron 
Department  of  Civil  Engineering 
Akron,  Ohio  44325 

Professor  Charles  W.  Bert 
University  of  Oklahoma 
School  of  Aerospace,  Mechanical, 
and  Nuclear  Engineering 
Norman,  Oklahoma  73019 

Professor  Satya  N.  Atluri 
Georgia  Institute  of  Technology 
School  of  Engineering  and 
Mechanics 

Atlanta,  Georgia  30332 

Professor  Graham  F.  Carey 
University  of  Texas  at  Austin 
Department  of  Aerospace  Engineering 
and  Engineering  Mechanics 
Austin,  Texas  78712 

Dr .  S .  S .  Wang 
University  of  Illinois 
Deoartment  of  Theoretical  and 
Applied  Mechanics 
Urbana,  Illinois  61801 

Professor  J.  F.  Abel 
Cornell  University 
Department  of  Theoretical 
and  Applied  Mechanics 
Ithaca,  New  York  14853 


Professor  V.  B.  Neubert 
Pennsylvania  State  University 
Department  of  Engineering  Science 
and  Mechanics 

University  Park,  Pennsylvania  16802 

Professor  A.  W.  Leissa 
Ohio  State  University 
Department  of  Engineering  Mechanics 
Columbus,  Ohio  43212 

Professor  C.  A.  Brebbia 
University  of  California,  Irvine 
Department  of  Civil  Engineering 
School  of  Engineering 
Irvine,  California  92717 

Dr.  George  T.  Hahn 
Vanderbilt  University 
Mechanical  Engineering  and 
Materials  Science 
Nashville,  Tennessee  37235 

Dean  Richard  H.  Gallagher 
University  of  Arizona 
College  of  Engineering 
Tucson,  Arizona  85721 

Professor  E.  F.  Rybicki 
The  University  of  Tulsa 
Department  of  Mechanical  Engineering 
Tulsa,  Oklahoma  74104 

Dr.  R.  Raftka 

Illinois  Institute  of  Technology 
Department  of  Mechanics  and  Mechanics 
and  Aerospace  Engineering 
Chicago,  Illinois  60616 

Professor  J.  G.  de  Oliveira 
Massachusetts  Institute  of  Technology 
Department  of  Ocean  Engineering 
77  Massachusetts  Avenue 
Cambridge,  Massachusetts  02139 

Dr.  Bernard  W.  Shaffer 
Polytechnic  Institute  of  New  York 
Route  110 

Farmingdale,  Nev  York  11735 
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Industry  and  Research  Institutes 


Industry  and  Research  Institutes  (Con’t) 


Dr.  Norman  Hobbs 
Kaaan  AviDvne 
Division  of  Raman 
Sciences  Corporation 
Burlington,  Massachusetts  01803 

Argonne  National  Laboratory 
Library  Services  Department 
9700  South  Cass  Avenue 
Argonne,  Illinois  60440 

Dr.  M.  C.  Junger 
Cambridge  Acoustical  Associates 
54  Rindge  Avenue  Extension 
Cambridge,  Massachusetts  02140 

Mr.  J.  H.  Torrance 
General  Dynamics  Corporation 
Electric  Boat  Division 
Groton,  Connecticut  06340 

Dr.  J.  E.  Greenspon 

J.  G.  Engineering  Research  Associates 
3831  Menlo  Drive 
Baltimore,  Maryland  21215 


Dr.  T.  L.  Geers 

Lockheed  Missiles  and  Space  Company 

3251  Hanover  Street 

Palo  Alto,  California  94304 

Mr.  William  Cayvood 
Applied  Physics  Laboratory 
Johns  Hopkins  Road 
Laurel,  Maryland  20810 

Dr.  Robert  E.  Dunham 
Pacifica  Technology 
P.0.  Box  148 

Del  Mar,  California  92014 

Dr.  M.  F.  Ranninen 
Battelle  Columbus  Laboratories 
505  Ring  Avenue 
Columbus,  Ohio  43201 

Dr.  A.  A.  Hochrein 
Daedalean  Associates,  Inc. 
Springlake  Research  Road 
15110  Frederick  Road 
Woodbine,  Maryland  21797 


Newport  News  Shipbuilding  and 
Dry  Dock  Company 
Library 

Newport  News,  Virginia  23607 


Dr.  James  W.  Jones 
Swanson  Service  Corporation 
P.0.  Box  5415 

Huntington  Beach,  California  92646 


Dr.  W.  F.  Borich 

McDonnell  Douglas  Corporation 

5301  Bolsa  Avenue 

Huntington  Beach,  California  92647 

Dr.  H.  N.  Abramson 
Southwest  Research  Institute 
8500  Culebra  Road 
San  Antonio,  Texas  78284 

Dr.  R.  C.  DeHart 
Southwest  Research  Institute 
8500  Culebra  Road 
San  Antonio,  Texas  78284 

Dr.  M.  L.  Baron 
Weidlinger  Associates 
110  East  59th  Street 
New  York,  New  York  10022 


Dr.  Robert  E.  Nickel 1 
Applied  Science  and  Technology 
3344  North  Torrey  Pines  Court 
Suite  220 

La  Jolla,  California  92037 

Dr.  Rev in  Thomas 
Westinghouse  Electric  Corp. 
Advanced  Reactors  Division 
P.  C.  Box  158 

Madison,  Pennsylvania  15663 

Dr.  B.  D.  Bibbitt 
Hibbitt  6  Rarlsson,  Inc. 

132  George  M.  Coban  Boulevard 
Providence,  Rhode  Island  02903 

Dr.  R.  D.  Mindlin 
89  Deer  Hill  Drive 
Ridgefield,  Connecticut  06877 
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Industry  and  Research  Institutes  (Con't) 


Dr.  Richard  E.  Dame 
Mega  Engineering 
11961  Tech  Koad 

Silver  Spring,  Maryland  20904 

Mr.  C.  M.  Stanley 
Lockheed  Palo  Alto  Research 
Laboratory 
3251  Hanover  Street 
Palo  Alto,  California  94304 

Mr.  R.  L.  Cloud 

Robert  L.  Cloud  Associates,  Inc. 
2972  Adeline  Street 
Berkeley,  California  94703 


V 


